Do the Unidentified EGRET Sources Trace Annihilating 
Dark Matter in the Local Group? 



J. Flix 1 *, J.E. Taylor 2 , M. Martinez 1 , F. Prada 3 , J. Silk 2 and J. 
Cortina 1 

Institut de Fisica d'Altes Energies, Universitat Autdnoma de Barcelona, E-08193 
Barcelona, Spain 

Astrophysics, University of Oxford, Denys Wilkinson Building, Keble Road, 
Oxford, 0X1 3RH, UK 

3 Instituto de Astrofisica de Andalucia (CSIC), E-18008 Granada, Spain 



Abstract. In a cold dark matter (CDM) framework of structure formation, the dark 
matter haloes around galaxies assemble through successive mergers with smaller 
haloes. This merging process is not completely efficient, and hundreds of surviving 
halo cores, or subhaloes, are expected to remain in orbit within the halo of a galaxy 
like the Milky Way. While the dozen visible satellites of the Milky Way may trace 
some of these subhaloes, the majority are currently undetected. A large number of 
high- velocity clouds (HVCs) of neutral hydrogen are observed around the Milky Way, 
and it is plausible that some of the HVCs may trace subhaloes undetected in the 
optical. Confirming the existence of concentrations of dark matter associated with 
even a few of the HVCs would represent a dramatic step forward in our attempts 
to understand the nature of dark matter. Supersymmetric (SUSY) extensions of 
the Standard Model of particle physics currently suggest neutralinos as a natural 
well-motivated candidate for the non-baryonic dark matter of the universe. If this 
is indeed the case, then it may be possible to detect dark matter indirectly as 
it annihilates into neutrinos, photons or positrons. In particular, the centres of 
subhaloes might show up as point sources in gamma-ray observations. In this work 
we consider the possibility that some of the unidentified EGRET 7-ray sources trace 
annihilating neutralino dark matter in the dark substructure of the Local Group. 
We compare the observed positions and fluxes of both the unidentified EGRET 
sources and the HVCs with the positions and fluxes predicted by a model of halo 
substructure, to determine to what extent any of these three populations could be 
associated. 
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1. Introduction 



Determining the nature of the unidentified EGRET sources (Hartman 
et al, 1999) is an important unsolved problem in modern high-energy 
Astrophysics. These sources constitute a major part of all the discrete 
sources sample detected by EGRET. Despite continuing efforts, most of 
them have not yet been associated with sources at other wavelengths. 
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Identification of EGRET sources is particularly difficult because of the 
large uncertainties in source positions and the lack of other observations 
at similar wavelengths. 

While a few of the unidentified EGRET sources may simply be arti- 
facts of the Galactic diffuse emission, it is normally thought that most 
are genuine detections, and that some fraction may even represent a 
completely new population of 7-ray emitters. As the sources are located 
predominantly at low to mid-Galactic latitudes, it seems likely they are 
mainly Galactic in origin, rather than being extragalactic objects like 
quasars (although this may account for some faction of them - e.g. 
(Sowards-Emmerd et al, 2003)). Thus most of the effort to identify the 
EGRET sources has focused on possible Galactic counterparts such 
as young pulsars (Hobbs et al, 2004), supernova remnants (Torres et 
al, 2003), isolated neutron stars (Goldoni and Mussol, 1996), molecu- 
lar clouds (Benaglia et al, 2001), objects in Gould's Belt (Gehrels et 
al, 2000), microquasars (Paredes, 2002), or young stellar objects like 
Wolf-Rayet stars or OB associations (Romero et al, 1999). 

On extragalactic scales, observational cosmology has recently con- 
verged on a 'concordance' model, which appears to satisfy all or most 
of the observational constraints simultaneously. In this 'ACDM' cos- 
mology, the present-day energy density of the universe is dominated by 
'dark energy', possibly in the form of a cosmological constant A, while 
the remaining matter content is constrained by nucleosynthesis and 
the microwave background to be predominantly non-baryonic cold dark 
matter (CDM) (Spergel et al, 2003). The nature of this dark matter 
is still unknown, and its identification is currently a top priority for 
theoretical and experimental Astrophysics. 

Weakly Interacting Massive Particles ( WIMPs) , non-relativistic relics 
from the early universe, have been proposed as one good candidate 
for non-baryonic dark matter. Although none of the particles in the 
Standard Model satisfy the requirements of being weakly interacting, 
massive and stable enough to be suitable relic particle, supersymmetric 
(SUSY) extensions of the Standard Model do provide a candidate that 
satisfies these requirements: the lightest SUSY particle (LSP), which 
ends up being one of the neutralinos (x) in most of the SUSY-breaking 
scenarios. The neutralino is a majorana particle and can annihilate with 
itself in pairs, producing high-energy neutrinos, positrons, antiprotons 
and 7-rays (Bergstrom et al, 2001). Thus in principle it can be detected 
through high-energy observations. 

In this paper we will discuss about the possibility that 7-ray emission 
from x self-annihilation, occurring in dark matter subhaloes in the halo 
of our Galaxy, may be responsible for some of the unidentified EGRET 
sources. We consider the spatial distribution and fluxes of these sources, 
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and compare them with the predictions of a model of halo substructure. 
We will also compare the source positions with the high-velocity clouds 
(HVCs), a population of neutral hydrogen clouds that may also trace 
dark matter concentrations in the galactic halo. 

2. Dark Matter Substructure in Galaxy Haloes 

In CDM cosmologies, galaxies form within dense concentrations of dark 
matter, called 'haloes' by analogy with the stellar halo of the Milky 
Way. Dark matter haloes are assembled through a process of accretion 
and the hierarchical merging of smaller subunits (cf. (Taylor et al, 
2003) and references therein). High-resolution numerical simulations 
of CDM structure formation have shown that this merging process is 
not completely efficient; they predict that hundreds of dense halo cores, 
or subhaloes, should remain in orbit around present-day galaxies (e.g. 
(Klypin et al, 1999)). While some of these subhaloes may correspond 
to the few dozen satellites which surround the Milky Way Galaxy and 
our neighbour M31 (Mateo, 1998), the expected number of satellites 
exceeds the number observed by approximately an order of magnitude 
(Klypin et al, 1999) (Moore et al, 1999). 

The current census of dwarf satellites is probably still incomplete. 
The local dwarf galaxies are often faint, low-brightness objects, and 
thus very difficult to detect, as evidenced by several discoveries made 
only recently (e.g. the Canis Major dwarf, (Martin et al, 2004)). Nonethe- 
less, it seems unlikely that most of the subhaloes predicted to orbit the 
Milky Way will ever be detectable by stellar emission alone. 

2.1. Subhaloes as 7-ray Sources 

The specific properties of the dark matter particle may provide an 
alternate means of detecting local substructure. Neutralino annihilation 
can generate continuum 7-ray emission, via the process \X ~^ QQ- The 
subsequent decay of 7r°-mesons created in the resulting quark jets pro- 
duces a continuum of 7-rays. Given a neutralino density profile pdm(r) 
(i.e. the dark matter density profile within a subhalo), the expected 
annihilation 7-ray flux above an energy Eq arriving to Earth is given 



by: 
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where iV 7 is the number of continuum gammas emitted per neu- 
tralino annihilation of m x mass above an energy Eq, and (crv) ann the 
averaged product of annihilation cross-section and relative velocities. 
The first part of Equation 1 can be evaluated for any given SUSY 
scenario, although the result varies by several orders of magnitude over 
the span of SUSY parameter space (Prada et al, 2004). 

All the information about the dark matter distribution is contained 
in the J(^) term. Cosmological simulations indicate that dark matter 
haloes are approximately spherical, with a cuspy central density profile 
that scales as Pdm( r ) °^ r~ a . The central slope of the density profile is 
thought to lie somewhere between a=l (Navarro et al, 2004) (the 'NFW 
profile') and a=1.5 (Moore et al, 1998), although recent simulations 
generally find that central slopes are shallower than 1.5. 

This profile should characterise both the main halo within which the 
Milky Way resides, and the subhalo cores that survive within it. The 
central cusp of the main halo may be affected by subsequent stages in 
galaxy formation, however, since baryons dominate the gravitational 
potential within the luminous part of bright galaxies. In particular, 
processes like starbursts, galaxy mergers or bar formation could heat 
and disrupt the dark matter cusp in the centre of the Milky Way, 
M31 or the brighter members of the Local Group. Smaller subhaloes 
should be less strongly affected by baryonic processes, and thus they 
may preserve 'pristine' dark matter cores, where the annihilation rate is 
still relatively high. Thus we will focus on comparing EGRET sources 
with small subhaloes, and ignore 7-ray emission (or lack thereof) from 
the brightest members of the Local Group. 

Using the semi-analytic model described in (Taylor et al, 2003), 
we have generated a number of realisations of halo substructure for a 
galaxy like the Milky Way. For each subhalo surviving at the present 
day, we calculated the 'flux multiplier' J(^) and the relative flux, as 
measured by an observer 8.5 kpc from the centre of the galaxy. A 
cut on subhalo masses of > 5 • 1O 7 M was applied, as the models are 
incomplete below this mass limit. 

In Figure 1 (left) we compare the predicted fluxes from the subhaloes 
in the each of our model systems with the fluxes measured for the 
unidentified EGRET sources. For each realisation, we have scaled the 
luminosity of the brightest subhalo to the luminosity of the brightest 
EGRET source. This gives an upper limit on the emissivity of the dark 
matter, for that particular model halo. By comparing the cumulative 
luminosity functions, we can also determine how many of the unidenti- 
fied sources could be subhaloes, assuming this scaling and a luminosity 
function like that of the model halo. This upper limit is between 10 
and 58 of the EGRET sources, depending on the realisation, with a 
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Figure 1. (left) Constrains to the maximum number of subhaloes compatible to 
unidentified EGRET sources, (right) Number distribution of all subhaloes with 
fluxes detectable by EGRET in a mSUGRA scenario. 



mean value of 26 ± 11. If we consider a neutralino of with m x = 100 
GeV, we can compute the number of gammas above the EGRET energy 
threshold with the Darksusy package (Gondolo et al, 2004). The (crv) gim 
needed to achieve such fluxes are from 1 • 10" 25 to 5 • 10" 24 cm s . 
These values exceed by one or two orders of magnitude the allowed 
WMAP (o"t>) ann in some SUSY scenarios. These numbers constrain an 
upper bound on the (<7t>) ann for a 100 GeV neutralino. 

As an example for a particular SUSY model, we consider to the 
set of mSUGRA scenarios discussed in (Prada et al, 2004). The range 
of models satisfying the WMAP relic density limits correspond to x 
masses between approximately 70 GeV and 1400 GeV, while {crv) gim 
lies in the range 1-10 -29 to 3T0~ 26 cm 3 s~ 1 . For these SUSY models, 
we compute the expected 7-ray flux above 100 GeV for all simulated 
subhaloes. The Figure 1 (right) shows the number distribution of all 
subhaloes with fluxes detectable by EGRET in this mSUGRA scenario. 
The number of detectable subhaloes vary from 1 to 5. Most favourable 
models correspond to those neutralinos with {crv) n ranging between 
3-10~ 26 cm 3 s _1 and 10 -27 cm 3 s _1 with masses from 100 to 400 GeV. 



3. Comparison with a Different Population of Halo Objects 

Given the small number of possible matches to between subhaloes and 
EGRET sources, it would be very useful to find another observable 
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tracer of halo substructure. The High- Velocity Clouds (HVCs) are iso- 
lated concentrations of neutral hydrogen whose velocities distinguish 
them from the gas in the galactic disk (Blitz et al, 1999). The kinematics 
of these clouds suggest that they may orbit the Milky Way at large 
distances, or that they may even be distributed throughout the Local 
Group. While the Compact HVCs (CHVCs) seem the most promising 
tracers of dark substructure (Braun and Burton, 1999), even the less 
compact clouds show similarities to subhaloes in their distribution, 
once cuts in halo distance are applied (Putman and Moore, 2001). The 
detection of a tidal interaction between a non-compact HVC and the 
dwarf galaxy LGS 3 further supports this association (Robishaw et al, 
2002). Thus we will compare the distribution of both types of clouds 
with that of our model subhaloes, and with the distribution of EGRET 
sources. 

3.1. An All-sky Catalogue of High Velocity Clouds 

Two different catalogues of High Velocity Clouds are available for the 
two hemispheres, and they can be combined to form an all-sky cata- 
logue of HVCs. The northern LDS survey (Heij et al, 2002) covered the 
sky north of declination -30°, while southern HIPASS survey (Putman 
et al, 2002) covered the sky south declination +2°. Both surveys are 
different in terms of sensitivity, angular resolution, velocity coverage 
and velocity resolution, but the same selection criteria was applied to 
both to identify isolated hydrogen clouds (Heij et al, 2002). In Table 3.1 
we list the main parameters of each catalogue. Approximate positional 
accuracies were estimated by comparing published values for known 
dwarf galaxies with the positions of the associated hydrogen clouds. 

Table I. Observational properties for HIPASS and LDS HVCs surveys. 



Vlsb, coverage 
V resolution 
5a T Sensitivity 
Angular resolution 
Positional accuracy 



LDS Survey 

-450 kms' 1 to +400 kms' 1 
1.0 kms- 1 
70 mK 
36' 
15' 



HIPASS Survey 

-700 kms' 1 to +500 kms' 1 
26.4 ferns" 1 
1 mK 
15.5' 
10' 



A total of 917 clouds were catalogued in the northern LDS survey, 
of which 777 are HVCs, 90 are CHVCs (isolated clouds in position and 
velocity), 37 were designated as CHVC: and 13 designated as CHVC? 
(clouds which were not unambiguously designated as CHVCs by the 
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Figure 2. Hydrogen clouds included in HIPASS survey. 
LPS HII Survey | 
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CHVC, CHVC: and CHVC? -» Zoom x 3.0 



Figure 3. Hydrogen clouds included in LDS survey. 

criteria described in (Putman et al, 2002)). For the HIPASS southern 
survey, a total of 1997 clouds were detected, of which 41 are actually 
galaxies (GLXY), 1618 are HVCs, 179 are CHVCs and 159 are desig- 
nated as :HVC as they could not be unambiguously classified as CHVCs 
because their elongation towards low level extended emissions. All the 
hydrogen clouds detected by either survey are shown in Figure 2 and 
Figure 3. 
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3.2. Positional Coincidences Between HVCs and 
Unidentified EGRET Sources 

Given the HVC positions from the HIPASS and LDS catalogues, we 
can test for positional coincidence between hydrogen clouds and 7- 
ray sources. More specifically, we can compare the positions of the 
maximum column density in neutral hydrogen for each cloud with the 
positions of the unidentified EGRET sources, on the supposition that 
both trace the dynamical centres of bound dark matter subhaloes. In 
order to determine the significance of positional coincidences, we define 
a relative distance between each hydrogen cloud and the each of the 
unidentified EGRET sources as: 

^2 l-^cgrct 2-cloudl 

°egrct + ^cloud 

where x egre t and x c ioud are the galactic longitude and latitude of the 
unidentified EGRET sources and the hydrogen clouds, expressed in 
degrees. We take as an characteristic uncertainty in the EGRET posi- 
tions degret = #95/2, where #95 is the 95% EGRET detection confidence 
contour. For cr c i ou d we choose the positional accuracies of HIPASS and 
LDS surveys given above. As <7 egr et > Ccioudi £ values will be dominated 

by Oegret- 




LDS HVC 



Figure 4- Unidentified EGRET sources whose positions coincide with HVCs at 
5 < 1, after applying 68% flux variability rejection. 

A total of 10 HIPASS HVCs and 2 LDS HVCs have an unidentified 
EGRET source lying within 5 <1 of their point of maximum emission. 
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As the number of unidentified EGRET sources in the field of view 
of HIPASS(LDS) surveys are quite similar, the different number of 
coincidences found for the two catalogues may be explained because the 
better sensitivity of HIPASS survey (HIPASS contains roughly double 
HI clouds than LDS) and worst LDS positional accuracies. 

The 7-ray emission from neutralino annihilation should be steady 
on observationally accessible timescales. The variability of all sources 
in the 3rd EGRET Catalogue has been evaluated, assuming a model 
of the flux distribution to characterise the flux variation, and using 
both positive flux detections and upper limits (Nolan et al, 2003). This 
method is quite successful at distinguishing between between the known 
(variable) quasars and the known (non-variable) pulsars. As we are only 
interested in those steady 7-ray emitters, independent of how much 
variable they are, we make use of the V12 parameter, which indicates 
the confidence with which the hypothesis of variability can be rejected 
for a given source (100% = steady source). If we apply a cut at the 
68% level to our 6 <1 sample, 6 sources are rejected, resulting in 5 
(1) remaining coincidences between HVCs to and unidentified EGRET 
sources for the HIPASS (LDS) surveys. These sources are shown in 
Figure 4 and listed in Table 3.2. 

To check for significance of these coincidences, we perform several 
simulations in which we distribute 10 4 fake EGRET sources isotrop- 
ically and evaluate the probability of finding coincident HVCs and 
EGRET sources, after including the rejection factor from variability. 
The average number of random associations found were 1.0 ± 0.3 for 
LDS survey and 3.1 ±0.5 for HIPASS survey, for the same coincidence 
cut criteria of 5 <1. Thus the real result obtained for HIPASS clouds is 
only marginally (4cr) above the level expected from random coincidence. 



Table II. Unidentified EGRET Sources and HVCs Found to be 
Spatially Coincident, f = HIPASS Survey, \ = LDS Survey. 



3rd EGRET Name 


HII Name 


a 


J0439+1105 


HVC 186.2-22.8-206 % 


0.18 


J0616-0720 


HVC 215.8-11.2+162 t 


0.54 


J1329-4602 


HVC 310.0+16.5+182 f 


0.51 


J1527-2358 


HVC 343.1+25.9-144 f 


0.98 


J1635-1751 


HVC 359.6+19.7-149 t 


0.32 


J1904-1124 


HVC 024.1-08.2-238 f 


0.48 
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4. Conclusions 

If the dark matter particle is a supersymmetric neutralino, one of 
the most popular candidates currently under consideration, then lo- 
cal concentrations of dark matter may be detectable as 7-ray sources 
distributed throughout the halo of the Milky Way. In this paper, we 
have compared the luminosities of unidentified sources in the EGRET 
7-ray catalogue with the luminosities of local dark matter substructure 
predicted by a semi-analytic model of neutralino annihilation in the 
halo of the Milky Way. We have also looked for spatial coincidences 
between EGRET sources and high-velocity neutral gas, which could be 
an observational tracer of local dark matter. 

We find that most of the unidentified EGRET sources cannot be 
dark matter subhaloes. Whatever the emissivity of annihilating dark 
matter, the cumulative luminosity function for the two populations is 
sufficiently different that only ~26 ± 11 of the EGRET sources could 
correspond to local substructure. Applying further cuts to sources on 
the basis of location or variability should place an even stronger con- 
straint on the number of possible sources, and thus on the properties 
of dark matter. 

We also find that most of the HVCs cannot correspond to EGRET 
sources either. While there are a few close positional coincidences be- 
tween the two populations, the overall number is roughly that ex- 
pected from random superposition. Thus there is no indication that the 
EGRET sources are associated with HVCs, whether or not the latter 
trace dark matter. This is interesting, because HVCs could also produce 
7-rays through other, less exotic mechanisms, such as interactions with 
cosmic rays. Comparing HVCs with sources in lower energy catalogues 
might produce a stronger constraint on this mechanism, however, since 
the spectrum of 7-rays produced in this case is broader than that 
expected for \ annihilation. 
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